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Abstract 

A new procedure of the linear position measurement which allows to obtain sensitivity 
better than the Standard Quantum Limit and close to the Energetic Quantum Limit is 
proposed and analyzed in detail. Proposed method is based on the principles of strobo- 
scopic quantum measurement and variation quantum measurement and allows to avoid 
main disadvantages of both these procedures. This method can be considered as a good 
candidate for use as a local position meter in the "intracavity" topologies of the laser 
gravitational-wave antennae. 

1 Introduction 

It has been known for many years that sensitivity of the weak force detectors based on "ordi- 
nary" position meters is limited by the Standard Quantum Limit Q (see discussion on "ordi- 
nary" and other position meters in [fj). On the other hand, the scientific society faced with the 
necessity to overcome this limit due to requirements of the gravitational- wave experiments ||. 
In order to detect gravitational force producing small relative displacement of the antenna's 
test masses, it is necessary to use methods that allow to overcome this fundamental limit of 
sensitivity. There are two promising measurement techniques that can be applied to solve this 
problem. These are stroboscopic measurement and variation measurement methods. 

Stroboscopic measurement J3J is the technique based on the fact that harmonic oscillator 
momentum perturbation at t = does not influence its coordinate at t + mr/uj m , where u m is 
the the oscillator's eigenfrequency and n is an integer. In other words, the back action noise 
does not influence the measurement result if instantaneous measurements are being made at 
regular time intervals of At = im/u m . Using the stroboscopic measurement it is possible to 
detect force 



where -Fsql is the value of force corresponding to the SQL and 6 is the duration of each 
measurement. 

Variation measurement technique can be divided into two main classes: the time-domain 
variation measurement || |(| and spectral variation measurement |7|]. Time-domain variation 
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measurement allows to eliminate the back action noise from the meter's output signal by in- 
troducing time-dependent cross-correlation between the meter's measurement and back-action 
noise. In the particular case of interferometric position meters being used in gravitational- wave 
detectors these noises correspond to phase and power shot noises of the optical pumping corre- 
spondingly, and this correlation can be introduced by modulation of the local oscillator phase. 
In the spectral variation measurement frequency- dependent cross-correlation is used. 

Unfortunately, all these methods have their own disadvantages. In order to achieve sensi- 
tivity substantially better than the SQL using interferometric position meter in stroboscopic 
regime it is necessary to use short pulses of optical pumping with high power W = irWo/uj m 9 3> 
Wo, where W is the mean power. Moreover, as 6 <C u;" 1 the meter's output signal contains 
additional noise due to brownian fluctuations of high-frequency internal modes of the test body. 
Estimates show that in the real gravitational-wave antennae this noise may be comparable or 
even larger than the fundamental quantum noise. Another disadvantage of the stroboscopic 
measurement is that it can be applied only on meters with the harmonic oscillator for the probe 
mass. At the same time only free masses are used for the test object in the contemporary laser 
gravitational wave antennae. 

On the other hand, in the case of the original form of time-domain variation measurement it 
is necessary to know shape and arrival time of the signal to be detected. As for the gravitational 
wave astronomy where the SQL overcoming is crucial, this procedure can not be applied because 
arrival time as well as exact shape of the signal are unknown. Spectral variation measurement 
is free from this disadvantage, but in the case of interferometric position meters it requires 
to use additional optical devices with bandwidth comparable with bandwidth of the main 
interferometer. For example, in the case of LIGO gravitational-wave observatory it is proposed 
to use additional Fabry-Perot cavities with the length comparable with the length of the main 
cavities of the antenna J71. 

In the articles |6], |8| the modified version of the time-domain variation measurement which al- 
lows to circumvent these disadvantages had been considered. This "discrete sampling variation 
measurement" (DSVM) is based on approximation of the signal by series of short rectangular 
"slices" with duration r < n/u max , where u mSuX is the upper signal frequency, and periodical 
applying of the variation measurement procedure fitted for such a rectangular pulse. In accor- 
dance with sampling theorem signal shape can be completely restored using signal values on 
consequent time intervals t n = nr as sampling coefficients. 

In this article we suggest new procedure which combines time-dependent pumping of the 
stroboscopic measurement, the meter noises time-dependent cross-correlation used in variation 
measurement, and the discrete sampling method. We propose to call this method "stroboscopic- 
variation measurement" (SVM). We will show that such a combined procedure allows to obtain 
sensitivity several times better than the pure DSVM procedure. This sensitivity can be close to 
traditional forms of the variation measurement and also close to the Energetic Quantum Limit 
0, [TO], which defines the ultimate sensitivity that can be achieved at given meter's energy. 

In the next section the idealized version of the SVM procedure based on the the instant 
position measurements is considered. In the section [5] we show that rigorous optimization of the 
variation measurements gives this procedure as a result. In the section ^ and in the appendix [A] 
we consider more realistic version of the SVM procedure based on position measurements with 
finite duration. In the appendix [FJ the Energetic Quantum Limit for linear position meters is 
considered in detail and it is also shown that sensitivity of the spectral variation measurement 
is equal to this limit. Appendix |C| devoted to the explicit solution of the optimization problem 
that arises in section [| 
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2 The idea of the SVM procedure 



Consider a slightly modified version of the triple measurement procedure discussed in the article 
IP- 1| - Let constant force F with duration r to act on a harmonic test oscillator with eigenfre- 
quency u m . In order to detect this force three instant measurements of the oscillator's position 
are performed at the moments t = 0, t = r/2, and t = r. The results of these measurements 
can be presented as 



Xj Xj -\- Xj fluct 



(2) 



where Xj (j = 1, 2, 3) are the "actual" values of the test object position at these moments of 
time, 
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mujr, 



mujr, 



(3a) 
(3b) 



where pi is the initial momentum of the test object and pi^ fluct are the test object momentum 
perturbations during the first and the second measurement correspondingly 

Estimate of the force F which does not depend on the initial state of the test object can be 
obtained using the formula 



muz 
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(4) 



Uncertainty of this value AF depends on the mean square errors for each of the measurement 
A X j, corresponding mean square perturbations of the test object momentum A p j, and the 
cross-correlations between the measurement errors and perturbations A xp j. They must satisfy 
the Heisenberg relation 
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(5) 



The experimentalist should optimize measurement by minimizing value of AF. Taking into 
account formula (|5|) it is easy to show that optimal values of cross-correlation are equal to 
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and in this case 
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Described procedure is very close to the variation measurement. Really, we eliminated the 
back-action term in the expression (0) using cross- correlation between 22fluct and p2fmct (see 
formulae @). In such a procedure it is possible to obtain, in principle, arbitrary high precision 
by increasing the values of A P j. 

On the other hand, if uj m T = irn then the term corresponding to the second measurement 
in the formulae (f|,|7|) vanishes. It means that only two measurements are necessary in this 
case, the first and the third ones, which correspond to the stroboscopic procedure. Therefore 
stroboscopic measurement which uses time- dependent pumping and variation measurement 
which uses time- dependent cross-correlation of noises should be considered as particular cases 
of "stroboscopic-variation" procedure described here. 

In the real experiment the values of A P j can not be arbitrary high due to energy limitations 
in the meter. Therefore value (|?|) can be further optimized taking into account additional 
condition that 



A 2 pl + A 2 p2 + A 2 p3 = S F T. 



(8) 



where Sf is some given value (sense of this peculiar notation will be explained below too). In 
the case of the interferometric position meter it is proportional to the optical energy stored in 
the interferometer, averaged over the time r. 

It is easy to show that in the optimal case we will obtain the following expression: 
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It is convenient to describe sensitivity by the meter effective noise spectral density 



as it was done in the article 
by the following formulae: 



S metcr = (AF) 2 r, (10) 
] . If the conditions (0) are valid, then the sensitivity is defined 
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is a convenient scale factor for the meter's noise equal to the Energetic Quantum Limit in the 
free test mass case taken at frequency (see Appendix ||). 

It should be noted that S^eter = *Seql only in the case of the pure stroboscopic measurement 
{uj m T = 7r), and S me ter > Seql for all other values of u m T (see the lowest curve in the Fig.[|). 

3 Linear position meter with time-dependent parame- 
ters 

Consider now more realistic situation when the meter which can continuously monitor position 
x(t) of the test object (for example, interferometric position meter) is used to detect signal 
force -FsignaiW- When the meter is turned on, its output signal can be represented as 



x(t) = x init {t) + £ fluct (t) + D-^FdgndCt) + F fluct (t)] , (14) 

where D is the linear differential operator describing the dynamics of the test object[| xn nc t(t) 
is the noise added by the meter, Fa nct (t) is the back-action force, and Xi n i t {t) is the "initial" 
position of the test object which describes its motion when signal force and back action force 
of the meter are absent. It should be noted that, whereas xa uct (t), Fa uct (t) and x init (t) are 
quantum-mechanical operators, x(t) is a classical observable (see discussion on de-quantization 



in linear position meters in the article [1 I 



Spectral densities of these noises S x and Sp, and their cross spectral density S x p satisfy the 
uncertainty relation 



2 ft 2 

S x Sp — S xF = — . (15) 

Suppose that the experimentalist can change values of these spectral densities within the bound- 
aries defined by the condition ([15]). If the interferometric position meter is used, the exper- 
imentalist can change the ratio Sf/S x changing the pumping power (as in the case of the 
stroboscopic measurement) and also he can change the value of S x f by changing the phase of 
the local oscillator (as in the case of the variational measurement). 

We suppose in this article for simplicity that noises Fa uct (t) and x init (t) can be regarded as 
"white" ones, and that the values of the spectral densities can be changed instantly (without 
any inertia). It means that the experimentalist have to be able to change the pumping energy 
in the interferometric position meter very quickly compared to the signal frequency. In the real 
experiments these approximations are valid only if relaxation time of the interferometric meter 
cavities is small compared to the signal characteristic time-scale (see brief discussion on this 
topic in the conclusion). 

On the first stage of signal processing information about the test object initial conditions 



(term x- m it(t)) is being excluded from the output signal x(t) by applying operator D to ( jl~4|) . 
As a result one will obtain the following expression for the estimation of the force F(t): 

In the case of the harmonic oscillator operator D can be written as 

d 2 , 
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F(t) = Dx(t) = F signal (t) + F fluct (t) + Dx fluct (t) 



(16) 



The next stage is the optimal estimation of the force F(t) averaged over short time interval, 
as it has been suggested in the DSVM method || ||. Mean-square error of this estimation is 
equal to 

t/2 t/2 

(AF)2 = T / { ~^§f dt+ J WtMt)I>v{t)+v(t)] 2 dt, (17) 

-t/2 -t/2 

where a = S x f/Sf is the variational factor and v(t) is the filter function which must satisfy 
conditions 

t/2 

J v{t) dt = l. (18) 

t/2 



t/2 

J S F (t)dt = S F r. (19) 

-t/2 

(compare with condition (§)). 

One can readily see that functional (fl7|) is a sum of two non-negative items, and variational 
coefficient a(t) appears only in the second one which originates in the back-action of the meter. 
It is evident that if variational factor a(t) satisfies the equation: 

a(t)Bv(t) + v(t) = (20) 
then the second term vanishes and expression ( |TTD takes form 



v(t) 



t=±r/2 



o, 



dv(t) 
dt 



t=±r/2 



o, 



We also suppose that function Spit) is limited by the condition 



-t/2 

Expressions (0), (|18|) , and (0) form the optimization problem to be solved in order to find 
characteristics for optimal signal processing. 

It should be noted that the standard Lagrange optimization procedure can not be used 
here, as function under integral in fl2"T| ) is not a continuously different iable one, so the Lagrange 
equations could not be solved. This problem can be solved using the Liapunov's problems 
optimization theory (See |L2| for details). The explicit solution is presented in Appendix 0, 



and here we will provide only the final expression. 
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(a) (b) 
Figure 1: Different layouts of the pumping pulses 



In the particular case of arbitrarily short pumping pulses the optimal function Sp{t) is equal 



to 



Sft 



s F (t) = {H 1 + ™^) 1 



6(t - r/2) + 2 cos ^f- S(t) + 6(t + r/2) 



if {jJ m T < 7T . 



if UJ m T > 71 , 



(22) 



The measurement error in this case is described by the formulae ( |TT| , P"2] ) . 

So we have shown that the sequence of the three instant measurements described in section 
^] represents an optimal procedure when the experimentalist is able to change freely the values 
of Sf, S x and S x f- 



4 Pumping pulses with finite duration 

It is evident, of course, that 5- like pumping pulses could not be obtained in practice. Therefore 
we consider here pumping pulses with finite duration. 

In order to reconstruct the signal force shape, the sequence of triple measurements described 
above have to be used. The experimentalist has a choice here, whether to use completely inde- 
pendent triads with each of them containing all three pulses (see Fig. 0(a)) or use overlapping 
triads with common first and third pulses (see Fig.|I](b)). The second variant is more "energy- 
saving" , but it does not permit to use variational technique for the first and the third pulses. 
It is possible to show that due to this reason its sensitivity is limited by the same condition ([!]) 
as sensitivity of the traditional stroboscopic procedure. Here we will consider the first variant 
only. 

In this case sensitivity is described by the formula 

_ h 2 m 2 uj _ (uj m r/7i) 4 

^meter — ~= . o~~/ \ — i 27 °EQL , 



where 



i r /2 

dj = = — / Sir(t) cos- 7 tu m t dt , where j = 1,2 . (24) 

Sf^ J-t/2 
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It depends on the shape of the pumping pulse, but it is evident that in any case the larger is 
Sf (i.e- the pumping power), the smaller is noise, and there are no absolute limitations similar 
to the SQL or the formula ([1]) here. 

Consider now quasi-optimal case when pumping pulses has small but finite duration 9 <C r. 
Suppose that 9 is the shortest pumping pulse duration that can be used in real experiment. In 
any case duration 9 should be large compared to the interferometer's cavities relaxation time. 
Let also the first and the third pulses to be divided between consequent measurement cycles, 
as it is shown in Fig. [I] (a). In this case function Sp(t) can be presented as 

S F {t) = S F r[kA(t + r/2) + (1 - fc)A(t) + kA(t - r/2)] , 

(25) 

where te [r (j - f) , r (j + §)] . 

Here factor k describes the relative energy of the pumping pulses, < k < 1, and A(t) is the 
narrow <5 like function that describes pumping pulse shape, and let 

r/2 r/2 r/2 

A(t)dt = l, J \t\ A(t) dt = a x 9 , J t 2 A(t)dt=(a 2 + a 2 1 )9 2 , (26) 

-r/2 -r/2 -r/2 

where aci, a 2 are some numeric factors which depend on the shape of the pulses. 

Substituting these expression into the formula ([^J), neglecting the terms of order of mag- 
nitude 9 3 or higher and optimizing the result with respect to k, we will obtain that 
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where T = r - 2a x 9 and f(9) w tt + (a\ - a 2 )uj 2 m 9 2 . 

These expressions correspond to the following optimal values of k: 



if u m T<f(9), 

if u m T>f(9), 
(27) 
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if u m T<f(6), 

if f(9)<u m T<ir 

if CUmT > 7T . 



(28) 



In the appendix [A] rectangular pumping pulses with arbitrary duration 9 are considered and 
explicit formulae for the measurement error are provided for this case. We suppose rectangular 
pulses to be a good model for real pumping pulses as the main disadvantage of square shaped 
function, that is infinite derivative on edges, does not influence the final expression for the 
measurement error, because due to fl2~4j) function Spit) that consists of rectangular pulses, is 
included to AF as integrand only and the final expression does not depend on its derivatives. 
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5 Conclusion 



Stroboscopic variation measurement presented in this article allows to obtain sensitivity better 
than the Standard Quantum Limit and close to the Energetic Quantum Limit using time- 
dependent values of meter noises spectral densities. In particular, if interferometric position 
meter is used then this procedure can be implemented by using both pumping power and phase 
of the local oscillator modulation. At the same time, it does not require pumping power to be 
in a non-classical state. 

It is important that the stroboscopic variation measurement does not require information 
about the signal shape and arrival time (as the usual variation measurement does). At the 
same time, sensitivity of the stroboscopic variation measurement does not depend crucially on 
the duration of the pumping pulses (as in the case of the usual stroboscopic measurement), see 

Fig. g. 

Evident area of this procedure application is the laser gravitational-wave antennae. How- 
ever, in authors' opinion, it can hardly be used in traditional topologies of gravitational-wave 
antennae due to two reasons. First, in this procedure the optical energy stored in the interfer- 
ometer must vary with characteristic time t <C Q^ 1 , where Q is the signal frequency. At the 
same time, relaxation time r* of the Fabry-Perot cavities used in the large-scale gravitational- 
wave antennae is close to the Q~ l which makes the energy modulation very difficult from the 
technological point of view. 

Second, in traditional topologies of gravitational- wave antennae very high values of the opti- 
cal pumping power are required in order to overcome the Standard Quantum Limit [ID] which, 



in particular, leads to undesirable non-linear effects in the large-scale Fabry-Perot cavities [13 



In the articles |Tg, [L6| a new class of so-called "intracavity" readout schemes for 
gravitational-wave antennae which allow in principle to increase sensitivity without increasing 
the pumping power is proposed. In such a schemes it is necessary to measure a small displace- 
ment of local test object with very high precision using, for example, a small (table-top scale) 
interferometric position meter. We think that stroboscopic variation measurement can be used 
most effectively in these meters. 
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A Rectangular pumping pulses 

Suppose pumping pulses we have introduced above to be square shaped. In this case function 
A(t) in ( P5D is equal to 

\ ' f 1(1 <S/2 ' (29) 
otherwise . 

Taking the above assumptions and formula (|T0| ) into account one can obtain the following 
expression for the measurement sensitivity: 
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where 



= 2 arccos - ( 2 cos 
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(32) 



Values 1 — k for different 0/t are plotted in Fig.|2]. 

Corresponding values of the measurement sensitivities are defined by the following formulae: 
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Figure 2: Graphics of 1 - k for 9/r =0, 0.2, 0.35, and 0.5. 



for the harmonic oscillator. 

These values are plotted in Figure |3|. Different curves correspond to several values of the 
pumping pulse duration 6. The topmost curve represents sensitivity of the DSVM procedure 
with constant pumping power ||. 

Here triangle ABC represents the "pure stroboscopic" area where k = 1. It should be noted 
that in the area right to this triangle the value of S^eter does not depend on pulse duration and 
remains the same as in the case of 5-function like pumping (|2l|). 



In order to detect the external force acting on some test object the following inequality 
should be fulfilled (see |TJ): 



where F(u) is the signal force spectrum, and S pos (u>) is the test object position fluctuations 
spectral density. 

Suppose position meter with back-action force spectral density Sf(uj) is attached to the test 
object. Let also this object dynamics to be described by function D(u), that is the operator 
D Fourier transform. Then spectral density S pos (uj) can be expressed in terms of back-action 
force spectral density Sf(uj) as 



B Energetic Quantum Limit 

This appendix is devoted to the formulae ([13]) derivation and also the expression for the ener- 
getic quantum limit when the harmonic oscillator is used as a probe body. 




(35) 



S f (oj) 



(36) 
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Therefore, condition ( p5|) can be rewritten as: 




(37) 
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Figure 3: Proposed method and DSVM procedure {£ = const) sensitivities (noise spectral 



densities) related to the EQL 



S*, 



motor 



EQL 



when 0/t =0, 0.2, 0.35, and 0.5. 



where 

ft 2 \D(co)\ 2 

W^) = T -^y- (38) 

Exactly the same condition can be obtained, if the uncertainty relation for linear position 
meter noises is used. Really, the total net noise of such a meter can be described by its spectral 
density 

StataiM = S F (u) + 2D{u)S xF {u) + \D(u)\ 2 S x (u) , (39) 
and spectral densities of the meter noises must satisfy the following condition: 

S x (u)S F (u) - S 2 xF (cv) > * . (40) 

Minimizing value Stotai(<^0 taking this additional condition into account, one can easily show 
that S totaX (u) > S eql (uj). 

Suppose that (i) Sp does not depend on frequency^ and (ii) we want to have S'eql 0^0 below 
some given value S'eql in some given spectral range < to < u; max . 

Let our test object be a free mass. In this case 

D(u) = -raw 2 , (41) 

and 

, . h 2 m 2 Ld 4 

SeqlH = ^-. (42) 
It is evident that this spectral density is maximal if u = tv max , and the maximum is equal to 

SeQL = SEQL(^max) = . G ■ (43) 

4op 



2 back-action noise is supposed to be "white" 
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This is the EQL for a free mass. 

Let the test object be a harmonic oscillator now. In this case we will obtain that 

D(u)=m(oj 2 m -u 2 ), (44) 

and 

h 2 m 2 (u 2 m -u 2 ) 2 

^EQL(^) = £gT • I 45 ) 

Suppose that we can freely tune u m in order to minimize the maximum value (s) of the Seql(^) 
in the spectral range < uj < u> ma , x . It is evident that in optimal case 

= % , (46) 
and maximum values of the Seql(w) are equal to 

SeQL = SeQl(O) = -SEQL^max) = max _ (47) 

loop 

This is the EQL for a harmonic oscillator. 

C The explicit optimization method 

In this appendix we will pay some attention to the explicit derivation of formulae presented 
above. 



Here we concentrate on the expression (|2lD method of optimization when conditions fllq) 
and ( |19D are applied. The harmonic oscillator case is considered, as the free mass case can be 
obtained simply by assuming uj m — ► 0. 

In order to reduce the two variables problem ([H]) into the single variable v(t) variational 
problem it is necessary to express function «Si?(t) in terms of v(t) using the Lagrange equation 
as follows: 

T^-A = o,=>*W = ^|iM«>l. (48) 

where A is the Lagrange multiplier. There is no difficulty now to transform the initial opti- 
mization problem into the Liapunov's one (See fl^ , |PTfl for details) by changing variable v(t) in 
accordance with the formula Dv(t) = u(t): 

Zr/2 I dt mill > Zr/2 M W dt = W m> ^gX 

J^ /2 sin(c; m (r/2 - £)) u(t) dt = 0, f^* cos(cj m (r/2 - t)) u(t) dt = 0. ' ' 



Using the ideas represented in [ 12" , |TT| the following minimum value of fl4"9|) can be obtained 



./^(t-|) + ^(t) + £5(*+f), w m re[0, tt] 
V^-^+T^ + f?) w ro r e (tt, oo), 
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where {/i, v, £} and {/3, 7, are indefinite coefficients that can be defined through substitution 
of function u(t) into the supplementary condition equations in ((49]) . The final expression for 
u(t) may be represented as 



- [5{t - I) - 2 cos ^<f(t) + 5(t + |)] , u m r e [0, tt] 
«(t) = < A f 0UB ^- (51) 
Y[*(*-^r) + *(*+5^)]. ^re(7r,oo). 

Filtering function t>(t) is expressed in terms of as 



1 



v(t) = — / sin(o; m (t - t'))u(t') dt! . 



2 



References 

V.B.Braginsky, Sov. Phys. JETP 26, 831 (1968). 



[2 

[3] 
[4] 
[5] 
[6] 
[7 



[9 
[10; 



[11 
[12 

[13 
[14 
[15 



V.B.Braginsky, M.L.Gorodetsky F.Ya.Khalili and K.S.Thorne, Physical Review D 61, 
4002 (2000). 

K.S.Thorne, 300 Years of Gravitation, Cambridge University Press, 1987. 
V.B.Braginsky, Yu.I.Vorontsov, F.Ya.Khalili, Sov. Phys. JETP-Lett. 33, 405 (1978). 
S. P. Vyatchanin E. A. Zubova, Physics Letters A 201, 269 (1995). 
S. P. Vyatchanin, Physics Letters A 239, 201 (1998). 

H.J.Kimble, Yu. Levin, A.B.Matsko, K.S.Thorne and S. P. Vyatchanin, Physical Review D 
65, 022002 (2002). 

S. L. Danilishin, F. Ya. Khalili and S. P. Vyatchanin, Physics Letters A 278, 123 (2000). 

V.B.Braginsky, F.Ya.Khalili, Quantum Measurement, Cambridge University Press, 1992. 

V.B.Braginsky, M.L.Gorodetsky, F.Ya.Khalili and K.S.Thorne, Energetic quantum limit 
in large-scale interferometers, in Gravitational waves. Third Edoardo Amaldi Conference, 
Pasadena, California 12-16 July, ed. S.Meshkov, Melville NY:AIP Conf. Proc. 523, pages 
180-189, 2000. 

V.B.Braginsky, M.L.Gorodetsky, F.Ya.Khalili, A.B. Matsko, K.S.Thorne, S. P. Vyatchanin, 
LANL preprint |gr-qc/ 01090031 submitted to Phys. Rev. D (2001). 



E. M. Galeev, V. M. Tikhomirov, The brief course of optimization problems, The MSU 
Publishing House, 1989. 

V.B.Braginsky S.E.Strigin and S. P. Vyatchanin, Physics Letters A 287, 331 (2001). 
V. B. Braginsky, F. Ya. Khalili, Physics Letters A 218, 167 (1996). 
V.B.Braginsky, M.L.Gorodetsky, F.Ya.Khalili, Physics Letters A 232, 340 (1997). 



14 



[16] V.B.Braginsky, M.L.Gorodetsky, F.Ya.Khalili, Physics Letters A 246, 485 (1998). 
[17] V. M. Alekseev, V. M. Tikhomirov, S. V. Fomin, The optimal control, Science, 1979. 



15 



